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ABSTRACT 



In wireless communications and navigation systems such as 
GPS, pseudorange measurements are made between a signal 
transmitter and a signal receiver based on pseudorandom 
noise codes modulating the transmitted signals. In the 
receiver, a code is tracked by forming a correlation function 
between the received code and a local code replica and 
locating the peak of the function using a delay-lock loop 
(DLL). Multipath signals, which reflect from objects before 
arriving at the receiver, degrade the code tracking perfor- 
mance. Multipath- invariant points are locations of the cor- 
relation function whose properties are independent of mul- 
tipath parameters, and they can therefore be located 
reproducibly in the presence of multipath. A multipath 
mitigation method consists of locating a multipath-invariant 
(MPI) point of an ideal autocorrelation function and mea- 
suring the distance between the MPI point and the DLL. The 
same MPI point is located in a received correlation function; 
and the distance between the point and the DLL, now 
affected by multipath, is measured. The difference between 
the ideal distance and the actual distance is the code tracking 
error resulting from multipath. The error is subtracted from 
the computed pseudorange or used to control the DLL. The 
method can be used to reduce the effects of all types of 
tracking error sources, such as signal transmission failure or 
code noise. 
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MULTIPATH AND TRACKING ERROR 
REDUCTION METHOD FOR SPREAD-SPECTRUM 
RECEIVERS 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

[0001] This application claims the benefit of U.S. Provi- 
sional Application No. 60/168,472 filed Dec. 1, 1999, which 
is herein incorporated by reference. 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 

[0002] This invention was supported in part by grant 
number 95-G-005 from the Federal Aviation Administration 
and the Department of Transportation (FAA/DOT). The U.S. 
Government has certain rights in the invention. 

FIELD OF THE INVENTION 

[0003] This invention relates generally to radio receivers. 
More particularly, it relates to techniques for reducing the 
effects of multipath in spread-spectrum radio communica- 
tions systems, such as the Global Positioning System (GPS). 

BACKGROUND ART 

[0004] The Global Positioning System (GPS) is a U.S. 
government-supported system for precise position and 
velocity determination of objects. GPS is used in both 
military and civilian applications such as aircraft navigation 
and surveying. Orbiting satellites emit coded radio fre- 
quency signals that are received and processed by receivers 
on or near the surface of the earth to obtain pseudorange 
measurements, approximate instantaneous distances 
between the satellites and receiver. Because the satellite 
orbits and signal transmission times are known precisely, the 
receipt time of a particular signal bit can be used to quantify 
the transit time or range to the particular satellite. The orbits 
of the GPS satellites are arranged in multiple planes so that 
signals can be received from at least four satellites at any 
point on or near the earth, allowing precise position and 
velocity measurements of the receiver. Each satellite con- 
tinually transmits two spread-spectrum signals that are 
modulated by at least one of two pseudo-random noise 
(PRN) codes unique to the satellite. A receiver can therefore 
identify and separate signals from each satellite into separate 
channels, process the channel data separately, and combine 
the processed data to compute the receiver position. In 
addition to the PRN code, the signals arc modulated by 
slower-varying data signals defining the satellite orbits and 
other relevant information needed for the computations. The 
two signal frequencies are in the L band and are referred to 
as LI and L2. The LI signal has a carrier frequency of 
1575.42 MHz, usually expressed as 1540f o , where f 0 =1.023 
MHz. The L2 signal has a carrier frequency of 1227.6 MHz, 
or 1200f o . 

[0005] The two PRN codes are known as the precise (P) 
code and the course/acquisition (C/A) code. The P code is 
available only to military and other authorized users, while 
the C/A code is available to all users for commercial and 
other civilian uses. The LI signal is modulated by both the 
P and C/A codes, and the L2 signal by the P code only. The 
P code is spread over a bandwidth of 20.46 MHz and has a 
repetition frequency of one cycle per week. The C/A code, 



which contains 1023 bits, has a 1.023 MHz clock rate and a 
repetition frequency of 1 kHz. The remainder of this descrip- 
tion focuses on the C/A code; however, analogous descrip- 
tions apply to the P code or any other PRN code. 

[0006] In the receiver, local signals corresponding to the 
known PRN codes are generated and correlated with the 
received signals. The phases of the local signals are adjusted 
until the local signals correctly track the received signal- 
Since the time at which each bit of the code sequence is 
transmitted from the satellite is known, the time of receipt of 
each bit is a direct measure of the transmit time of the signal 
from the satellite to the receiver, and therefore a measure of 
the distance between the two. Based on the computed 
relative phase of the received signal, the receiver calculates 
the desired quantities of distance, velocity, etc. 

[0007] In order to understand the context of the present 
invention, it is necessary to understand the method by which 
the receiver tracks the incoming code. The most common 
pseudorandom noise code is a binary sequence of 0's and 1 's 
or -l's and +l's that modulates the carrier phase. The code 
sequence appears to be random, but in fact repeats identi- 
cally at regular intervals. The C/A code consists of 1023 
"chips," the length of time during which the code remains 
constant; the C/A chip period is therefore approximately 
977.5 nanoseconds, equivalent to a distance of 293.05 
meters. A portion 10 of a C/A code is shown in FIG. 1. In 
the receiver, a correlation function is formed between the 
actual received signal (after filtering) and a local replica. 
FIG. 1 illustrates the concept of the correlation function 
with code replicas 12a-12/ and 14 that are offset from 
received code portion 10 by whole numbers of chips. Each 
chip period in which the received and replica value of the 
code are identical contributes +¥n (the number of chips in 
this example) to the correlation function, and each chip with 
different values of the codes contributes -V12. When the 
replica is directly lined up with the code, as is the final code 
replica 14 of FIG. 1, the correlation function has a value of 
1. In all other offsets, as for replicas 12a-12f, the correlation 
function has a magnitude much less than 1. For longer 
sequences, the correlation function is much smaller in magr 
nitude and is often considered to be effectively zero beyond 
an offset of one chip in either direction. 

[0008] The bulk of the power in the correlation function 
occurs within a region between +1 and -1 chips offset from 
the received code signal. FIG. 2 shows a received code 
signal 16 and two offset replicas 18 and 20, along with a 
correlation function 22. The early 26 and late 28 arrows on 
function 22 correspond to code replicas 18 and 20, respec- 
tively. Clearly, offsetting replicas 18 and 20 from the 
received function by a value less than one chip decreases 
correlation function 22 from its maximum value of one, but 
still provides enough overlap for it to be significantly greater 
than zero. As a result, correlation function 22 has a triangular 
central portion 24 referred to as its main lobe. A standard 
delay-lock loop (DLL) uses a pair of correlators 26 and 28 
that sample function 22 at early and late offsets from the 
peak 30. The goal of the DLL is to straddle the peak 
symmetrically so that peak 30 is exactly half-way between 
correlators 26 and 28. Practically, correlation function 22 is 
sampled at the early and late positions 26 and 28, and an 
early minus late value is obtained and used to control the 
offsets of the correlators. Peak 30 is considered to be located 
when the early minus late value is zero. If the early minus 
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late value is positive, then the correlators are shifted earlier. 
Analogously, if the early minus late value is negative, then 
the correlators are shifted later. 

[0009] A highly schematic view of a GPS receiver 31 is 
shown in FIG. 3. The RF front end of most conventional 
GPS receivers is analog. An incoming signal enters through 
an antenna 32 and is immediately amplified by a pre- 
amplifier. The signal is then cable-routed to a multi-stage 
downconversion process 34, during which a series of inter- 
mediate frequency mixers and filters converts the signal, 
nominally 20 MHz wide and centered at 1.545 GHz for LI, 
C/A code receivers, to baseband for correlation, tracking, 
and subsequent processing. A final intermediate frequency 
stage filter 36 sets the precorrelation bandwidth (PCB). 
Narrowband receivers have a PCB less than 2.5 MHz wide 
and pass only the main lobe of the C/A code power spec- 
trum. Some implications of narrow PCB receivers are dis- 
cussed below. The signal is then digitally sampled (38), and 
the values of the correlation function at the late (40) and 
early (42) correlators are determined. A processor 44 com- 
putes the early minus late error and determines the required 
adjustment to the correlators. A digitally-controlled (or 
numerically-controlled) oscillator (DCO) 46 adjusts the 
phase of the early 48 and late 50 code replicas until the code 
is correctly tracked. 

[0010] This standard DLL implementation assumes a 
direct line of sight (LOS) between each satellite and the GPS 
receiver, i.e., that the relative time delay between signal 
transmission and reception is equal to the time for the signal 
to travel directly from the satellite to the receiver. In fact, 
most received signals are distorted by reflections of the 
transmitted signal from objects such as buildings and moun- 
tains or from the ground. Thus the received signal is a sum 
of the LOS signal and the so-called multipath (MP) signals. 
A highly schematic illustration of multipath is shown in 
FIG. 4. A signal transmitter 52, e.g., a satellite, transmits a 
signal that is received by a receiver 54. A line of sight signal 
56 travels directly from transmitter 52 to receiver 54. 
Additionally, multipath signals 58 and 60 reflect from a 
mountain 62 and building 64, respectively, before arriving at 
receiver 54. While only two multipath signals are shown, 
there are typically many more multipath signals arriving at 
receiver 54. Multipath signals are always delayed relative to 
the LOS signal, because they necessarily travel a longer 
distance between the transmitter and receiver. As a result of 
the reflections, multipath signals typically also have lower 
amplitudes than LOS signals. 

[0011] The effect of multipath on the DLL is illustrated 
schematically in FIG. 5. For clarity, only one multipath 
signal is shown. A multipath code 66 (dashed fine) adds to 
a line-of-sight signal 68 (dotted line) to yield a total code 
signal 70 (solid line). Note that MP signal 66 is a delayed 
and lower- amplitude version of LOS signal 68 and therefore 
adds a delayed, lower-amplitude function 72 to the LOS 
correlation function 74. As a result, the peak shape of the 
total correlation function 76 is distorted, causing the corr- 
elators to converge to a location 78 that is offset from the line 
of sight peak 80. The difference between LOS peak 80 and 
the midpoint between the two correlators (often tracked with 
an additional "prompt" correlator) is the DLL tracking error 
82 caused by multipath. The tracking error in the pseudor- 



ange measurement propagates to position, velocity, and 
other measurements, and is therefore highly detrimental to 
the accuracy of the system. 

[0012] Multipath tracking errors are difficult to remove for 
a number of reasons. Multipath errors are not zero mean* 
particularly for large amplitude MP signals, so that even 
infinite smoothing of the computed pseudorange cannot 
guarantee unbiased position errors. In general, the pseudo- 
range is a nonlinear function of multipath parameters such as 
amplitude, delay, and phase, which are not constant; changes 
in any one of these parameters can affect the DLL in 
unpredictable ways. Multipath errors are also particular to 
the location of the receiver, and thus errors at one receiver 
cannot be used to eliminate errors at a second receiver. 

[0013] The multipath tracking error can be reduced some- 
what by decreasing the spacing between the early and late 
correlators. Multipath effects occur primarily at the late 
slope of the main lobe, particularly for long-delay multipath. 
Thus while traditional DLL early and late correlators have 
spacings of one chip, many current receivers have correlator 
spacings of 0.5 chips or less. Decreasing the correlator 
spacing, however, cannot eliminate the fundamental effects 
of multipath. Furthermore, receivers with a narrow PCB 
flatten the correlation peak, and narrow correlator spacing 
has little effect on the tracking error in such receivers. FIG; 
6 is a plot of theoretical tracking error as a function of 
multipath delay for two different correlator spacings, 0.1 and 
1 chip, for a signal-to-multipath ratio of 3 dB and a 2 MHz 
precorrelation bandwidth. 

[0014] A large number of techniques for eliminating or 
reducing multipath-induced tracking error have been devel- 
oped. Current techniques can be grouped into two major 
categories, separation and estimation. Separation techniques 
essentially attempt to separate the LOS and MP signals and 
then track the LOS signal alone to eliminate the effects of 
multipath. Estimation techniques estimate the parameters of 
the LOS and MP signals to approximate their combined 
effect on the tracking error. All of these techniques require 
extra hardware for additional processing of the correlation 
function and often a wider bandwidth. Typically, additional 
correlators sample the peak at many locations to map out the 
complete curve. These methods are sensitive to changes in 
multipath parameters and have difficulty as the multipath 
parameters approach the line of sight parameters, i.e., at low 
multipath delays. Under these conditions, distinguishing 
LOS from multipath requires a very precise sampling of the 
function. Most multipath mitigation techniques therefore 
perform well against long-delay multipath, but for delays of 
below 30 meters, their performance degrades to that of a 
conventional, one-chip spacing DLL. Current techniques are 
also primarily effective only in wideband receivers. Receiv- 
ers with narrow PCB generate correlation functions with 
flattened and rounded main lobes; even with relatively long 
multipath delays, distinguishing multipath from LOS is 
difficult in such receivers. 

[0015] The E1-E2 Tracker technique is disclosed in P. G: 
Mattos, "Multipath Elimination for the Low-Cost Consumer 
GPS "Proceedings of the 9 th International Technical Meet- 
ing of the Satellite Division of the Institute of Navigation, 
Vol. 1, pp. 665-671, 1996. In this method, an additional pair 
of correlators samples the main lobe of the correlation 
function at its earliest, relatively undelayed portions. The 
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computed actual slope obtained from the correlators is 
compared with an ideal line of sight slope, and the difference 
is used to control the DLL correlators dynamically. This 
technique has a number of drawbacks. For accurate esti- 
mates of the slope, a maximum receiver precorrelation 
bandwidth is required. That is, it does not work well in 
narrow PCB receivers. It also performs best against 
medium- and long-delay multipath, and performs poorly 
under high noise conditions. Thus the E1-E2 Tracker method 
has similar drawbacks to other multipath mitigation tech- 
niques. 

[0016] It would be a significant improvement in the art to 
provide a multipath mitigation technique that is effective in 
narrow PCB receivers. Narrowband receivers provide a 
number of advantages over wideband receivers. They are 
much less vulnerable to both narrowband interference and 
degradations or anomalies in the transmitted signals. It 
would also be an improvement to provide a multipath 
mitigation technique that operates independently of multi- 
path parameters. 

OBJECTS AND ADVANTAGES 

[0017] Accordingly, it is a primary object of the present 
invention to provide a multipath mitigation technique that: 

[0018] 1) performs well in narrow precorrelation band- 
width receivers; 

[0019] 2) performs well with short -delay multipath; 

[0020] 3) operates independently of multipath param- 
eters; 

[0021] 4) requires few changes to existing receiver 
hardware; and 

[0022] 5) can be implemented with very low additional 
cost. 

[0023] It is an additional object of the present invention to 
provide a method for measuring a code tracking error 
generated by a spread-spectrum receiver tracking a coded 
signal. 

SUMMARY 

[0024] These objects and advantages are attained by a 
method for mitigating the effect of multipath signals on code 
or envelope tracking of a primary (line-of-sight) signal by a 
spread-spectrum receiver such as a Global Positioning Sys- 
tem (GPS) receiver. The primary signal, typically the earliest 
arriving signal, is a carrier signal modulated by an envelope, 
and the receiver tracks the envelope. The receiver receives 
a signal that consists of the primary signal plus the multipath 
signals and attempts to track the signal in order to compute 
a pseudorange measurement between the receiver and a 
signal transmitter. The method has the following steps: in an 
ideal comparison (autocorrelation) function, locating an 
ideal multipath-invariant point, thereby determining an ideal 
distance between the ideal multipath-invariant point and an 
ideal multipath-dependent point; in an actual comparison 
(correlation) function computed from the received signal, 
locating a corresponding multipath-invariant point, thereby 
determining an actual distance between the corresponding 
multipath-invariant point and a corresponding multipath- 
dependent point; comparing the ideal and actual distances; 
and using this comparison to mitigate the effect of the 



multipath signals on envelope or code tracking of the 
primary signal. In one embodiment, the difference between 
the actual and ideal distances (the tracking error) is sub-^ 
traded from the pseudorange measurement. In an alternative 
embodiment, the tracking error is used to control a signal 
tracking or delay-lock loop that tracks the received signal; 
The receiver typically also tracks a number of additional 
signals, and the tracking errors of all signals can be averaged 
to yield an average tracking error. The method preferably 
also includes a step of removing a bias in the tracking error; 

[0025] Preferably, the ideal multipath-invariant point is in 
a plateau of the ideal autocorrelation function, and most 
preferably at an edge of the plateau. Practically, it is useful 
to select the multipath-invariant point at a location that is 
adjacent to the plateau but not within it. The multipath- 
dependent point is typically at one of the sampling locations 
of the signal tracking loop. 

[0026] A multipath-invariant point in the actual correla- 
tion function is preferably found by bounding an interval in 
which the point is likely to occur, based on theoretical 
tracking errors and the ideal distance, and narrowing the 
interval until the point is located to a sufficient precision. If 
only a few samples (e.g., correlator observation or measure- 
ment points) of the multipath-invariant region are available, 
the interval is preferably narrowed according to a modified 
Fibonacci search. Preferably, the correlation function is 
sampled at a minimum of two locations, and a difference is 
formed between the two samples to obtain a discriminator 
value that is compared with a predetermined threshold 
value. Sample locations are adjusted until the discriminator 
value matches the threshold. Alternatively, a curve fit of the 
function is generated from the sampled points, and the 
multipath-invariant point is the point yielding a predeter-l 
mined threshold value of the curve. The method can also 
locate more than one multipath invariant point, and the 
tracking errors obtained from all points can be averaged to 
obtain an average tracking error. Preferably, the effect of 
signal tracking loop dynamics is removed from the raultir 
path sampling use a low-distortion component of the signal 
such as a carrier signal. Also, the correlation function is 
preferably normalized or calibrated before the multipart^ 
invariant point is located. 

[0027] The present invention also provides a method for 
measuring tracking error generated by a receiver that tracks 
a received signal modulated by an envelope. This method 
has steps similar to the multipath-mitigation method, but 
instead locates points that are invariant to signal distortion, 
The difference between the actual and ideal distances is the 
tracking error. For example, the tracking error can be caused 
by failure in signal transmission or by noise in the received 
signal. The measured tracking error can be used to improve 
performance of the receiver or to detect the presence of such 
failures of the signal transmission. The method can also 
include removing a bias in the tracking error, calibrating or 
normalizing the correlation function, or decoupling sample 
dynamics and noise from the distortion-invariant point loca4 
tion process using a low-distortion component of the signal 
such as a carrier signal. 
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BRIEF DESCRIPTION OF THE FIGURES 

[0028] FIG. 1 (prior art) illustrates the computation of an 
autocorrelation function of a small code segment. 

[0029] FIG. 2 (prior art) shows an ideal autocorrelation 
function with a pair of delay-lock loop correlators. 

[0030] FIG. 3 (prior art) is a block diagram of part of a 
conventional GPS receiver. 

[0031] FIG. 4 (prior art) is a schematic illustration of 
line-of-sight and multipath signals arriving at a receiver. 

[0032] FIG. 5 (prior art) illustrates the effect of multipath 
on a correlation function of a received signal and on code 
tracking by a delay-lock loop. 

[0033] FIG. 6 (prior art) is a plot of tracking error versus 
multipath delay for a conventional narrow precorrelation 
bandwidth GPS receiver, for two different correlator spac- 
ings. 

[0034] FIGS. 7A and 7B are ideal autocorrelation func- 
tions for PRN 8 and 29, respectively. 

[0035] FIG. 8 shows multipath-invariant regions of an 
ideal autocorrelation function. 

[0036] FIG. 9 illustrates a method of the present invention 
for eliminating multipath tracking error. 

[0037] FIG. 10 is a flow diagram of a method of the 
present invention for eliminating multipath tracking error. 

[0038] FIG. 11 illustrates the error introduced by narrow 
precorrelation bandwidth receivers in locating a multipath- 
invariant point. 

[0039] FIG. 12 is a block diagram of part of a receiver for 
implementing the present invention. 

[0040] FIG. 13 is a block diagram of the Tracking Error 
Compensator (TrEC) of the receiver of FIG. 12. 

[0041] FIGS. 14A-14C illustrate details of one implemen- 
tation of the method of FIG. 10. 

[0042] FIG. 15 illustrates a multiple correlator curve fit 
embodiment for locating a multipath-invariant point. 

[0043] FIG. 16 illustrates a multiple correlator embodi- 
ment for locating multiple multipath-invariant points. 

[0044] FIG. 17 is a plot of tracking error versus multipath 
delay for conventional narrowband and wideband multipath 
mitigation techniques and for a method of the invention. 

DETAILED DESCRIPTION 

[0045] Although the following detailed description con- 
tains many specifics for the purposes of illustration, anyone 
of ordinary skill in the art will appreciate that many varia- 
tions and alterations to the following details are within the 
scope of the invention. Accordingly, the following embodi- 
ments of the invention are set forth without any loss of 
generality to, and without imposing limitations upon, the 
claimed invention. 

[0046] The present invention provides a novel technique 
for mitigating the effects of multipath signals on code 
tracking by spread -spectrum receivers. While the invention 
will be described with respect to Global Positioning System 
(GPS) receivers, it is to be understood that the principles can 



be applied to any wireless information transfer or commu- 
nications systems (e.g., spread spectrum radio communica- 
tions systems) in which tracking of a carrier signal whose 
amplitude or phase is modulated by an envelope is per- 
formed. In addition, while multipath mitigation is a key 
application of the invention, other applications benefit from 
the ability of the invention to discern the difference between 
an ideal, uncorrupted signal and a signal distorted by mulr 
tipath, signal generator failure, or other sources. That is, the 
present invention can be used to identify and eliminate 
tracking errors originating from virtually any source. 

[0047] In the following description, a number of specific 
terms are used to illustrate the invention. It is to be under- 
stood that these terms are purely illustrative, and that inven- 
tion is not limited to such specific terms. The term "con> 
elator" or "correlator pair" describes a receiver component 
responsible for sampling a correlation function at particular 
locations. These terms apply to conventional analog receiv- 
ers. In all -digital (i.e., direct-sampling) receivers, however, 
the same function is performed by sampling the digital 
correlation function at chosen locations. It is to be under- 
stood that the use of the term "correlator" below applies both 
to analog correlators and to sampling the correlation funcr 
tion in digital receivers, and that its use does not limit the 
scope of the invention to any particular type of receiver. : 

[0048] A code is a specific example of an envelope that 
modulates the amplitude or phase of a carrier signal. In the 
context of the present invention, a receiver is tracking the 
envelope, and the invention detects and eliminates errors in 
envelope tracking. This process is typically referred to as 
code tracking. A delay- lock loop (DLL) is a specific receiver 
component that acquires and tracks an incoming signal. In 
the context of the present invention, the DLL tracks the 
envelope modulating a carrier signal. This function can be 
performed by any wireless communication or information 
transfer signal tracking loop, referred to herein as a "signal 
tracking loop." Any feedback signal control loop can be 
used, and it is to be understood that the term DLL does not 
limit the scope of the invention. As used herein, a primary 
signal refers to a signal that is tracked by a receiver in the 
presence of multipath or other distortions that cause tracking 
error. Typically, it is the earliest arriving signal at a receiver; 
The primary signal is referred to below as the line-of-sight 
(LOS) signal, but it is to be understood that the present 
invention improves tracking of any type of primary signal: 

[0049] As used herein, the term "comparison function" 
refers to any function in which a specific type of multipLU 
cation is performed between two signals. In particular, at 
each possible offset of the two signals, corresponding 
samples of the two signals are multiplied, and then all such 
products are summed to obtain a comparison function value 
for that offset. In GPS receivers, the comparison function is 
referred to as a correlation function that is generated 
between a received signal and, for example, a local replica 
of the received signal. The term comparison function 
includes, but is not limited to, a matched filter operation; 
convolution, integrate- and -dump accumulator, and envelope 
detector. 

[0050] As described above, it is generally assumed that the 
value of a pseudo-random noise (PRN) code correlation 
function is effectively zero outside of the main peak. Howt 
ever, the PRN codes are not strictly orthogonal. That is, the 
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correlation function is not identically zero away from the 
main lobe. The present invention arises from the insight that 
side lobes of the correlation function have valuable infor- 
mation that can be used to verify the accuracy of the 
delay-lock loop. While the main lobe is almost always 
affected by multipath, there exist regions of the correlation 
function, which may lie away from the main lobe, whose 
properties are not affected by multipath. If these regions can 
be located reproducibly, regardless of multipath effects, then 
they can be used to correct the tracking error in the delay- 
lock loop caused by multipath. It is commonly believed that 
the signal power outside of the main lobe is much too low 
to provide useful information. The present invention, how- 
ever, has been shown to extract useful information from the 
side lobes and use this information to eliminate or reduce 
multipath-induced tracking error. 

[0051] Locations of the correlation function whose prop- 
erties, when sampled and averaged over a sufficient length of 
time, are independent of multipath parameters such as 
amplitude, phase, and delay are called multipath-invariant 
(MPl) regions or points. Each PRN code has a unique set of 
multipath-invariant points. One useful multipath-invariant 
property is the slope of the correlation function. For illus- 
tration purposes, the invention is described below in the 
context of this property. 

[0052] Multipath invariance is best understood by consid- 
ering the two autocorrelation functions shown in FIGS. 7A 
and 7B for PRN 8 and PRN 29, respectively. (Recall that 
each GPS satellite transmits signals modulated by a unique 
and uniquely-numbered PRN code.) The autocorrelation 
function for PRN 8 contains side lobes 84a and 84b that are 
adjacent to a main lobe 86. Delay-lock loops (DLLs) track- 
ing PRN 8 have one of the worst multipath performances, 
because side lobe 84a overlaps main peak 86 even for 
relative multipath delays of as high as 3.5 chips. Conven- 
tionally, it is assumed that the tracking error is zero for 
relative multipath delays greater than 1.5 chips, because the 
main lobe of the multipath signal no longer interferes with 
the main lobe of the line of sight signal. In the case of PRN 
8, however, side lobe 84a interferes with DLL tracking just 
as a small replica of the main peak would, even if the main 
peak itself does not interfere. PRN 29, in contrast, has long 
plateaus 92a and 92b immediately adjacent to the main lobe 
94. In this case, multipath delays of greater than 1.5 chips 
have virtually no effect on DLL tracking, regardless of the 
amplitude, phase, and phase rate characteristics of the mul- 
tipath. 

[0053] in fact, if it were possible to track a side lobe 88 of 
PRN 8 located at -5 chips (5 chips before the main peak), 
side lobe 88 would display the multipath tracking error 
characteristic of peak 94 of PRN 29, i.e., virtually no 
tracking error. Multipath signals with relative delays of 
greater than 1.5 chips but less than the distance between side 
lobe 88 and its adjacent side lobe 90 have no effect on the 
slope of side lobe 88 (although minor changes in amplitude 
do occur). In general, regions of a correlation function that 
are not affected by multipath delays of up to some large 
number of chips can be considered to be multipath invariant. 

[0054] Furthermore, locating these regions, at least theo- 
retically, is relatively straightforward. FIG. 8 shows a 
generic correlation function with multipath-invariant 
regions 96a-96e indicated. Multipath regions 96a-96e are 



plateaus in the correlation function, i.e., regions with essen- 
tially zero slope. Superimposing a zero-slope multipath 
component onto a zero-slope line-of-sight component 
changes the amplitude, but not the slope, of the overall 
signal. In particular, regions 96a-96e are at the early part of 
the plateaus. For example, region 96b begins (i.e., on the 
right) at a point 98 with a discontinuous slope (between zero 
and non-zero slope regions) and extends until a point 100 a 
fixed distance 102 away from a later point 104 of discon- 
tinuous slope. Fixed distance 102 can be selected based on 
an expected or reasonable multipath delay. The selection of 
distance 102 assumes that the multipath signal will not be 
delayed so much that the lobe beginning with point 104 
interferes with region 96b. 

[0055] Practically, it is necessary not just to locate a 
multipath-invariant region, but to locate one (or more) 
multipath-invariant point. For example, the simplest such 
point to locate that is within MPI region 96b is point 98. In 
general, reproducibly locatable MPI points occur at the 
corners or slope transition points alongside the early slopes 
of the side and main lobes, i.e., at the edges of the plateaus; 
Multipath-invariant points are selected by generating an 
ideal autocorrelation function for each PRN code and iden- 
tifying plateaus within the function that are longer than a 
selected value. For example, it is reasonable to assume that 
multipath having relative delays (i.e., relative to the LOS 
signal, not relative to the satellite transmission) of at least 
three chips, approximately 880 meters, are either highly 
unlikely or have extremely low mean signal-to-noise ratios, 
compared with signal-to-noise ratios of the LOS main and 
side lobes. This assumption is particularly valid when the 
receiver and reflector geometry is changing, i.e., when the 
receiver or reflector is moving, but is also valid for large 
distances between the reflector and receiver, when higher 
frequencies have a more substantial decrease in amplitude '■. 
Both factors tend to "whiten" the MP interference, bringing 
it closer to zero mean and therefore easier to filter. 

[0056] Table 1 lists optimal multipath-invariant point 
locations, relative to the peak of the main lobe, and widths 
of the corresponding plateau for each GPS Gold Code PRN; 
Only 17 of the codes have multipath-invariant plateaus 
adjacent to the main lobe. Most multipath-invariant plateaus 
have a normalized plateau height of - 1 /i023, but they also 
occur atop trapezoidal side lobes. The minimum width of the 
optimal MPI plateau for all codes is ten chips wide (for PRN 
20). 
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TABLE 1-continued 
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[0057] A multipath -invariant point is used as shown in 
FIG. 9 to mitigate the effects of multipath on the line-of- 
sight code tracking performance of the delay-lock loop of 
the receiver. An ideal distance 106 between a selected 
multipath-invariant point 108 and, for example, the early 
correlator 110 of the delay-lock loop is measured on an ideal 
autocorrelation function generated from the PRN code. Ideal 
distance 106 is the expected distance for a strictly line-of- 
sight signal, i.e., one that does not include multipath signals 
(or other significant error sources). In the receiver, an 
additional correlator pair locates the corresponding multi- 
path-invariant point in the received correlation function. 
Within a reasonable error, which depends on receiver band- 
width, the corresponding point in the received signal is 
assumed to be at the same location, relative to the LOS peak, 
as the ideal multipath-invariant point in the ideal autocor- 
relation function. An actual or measured distance 112 
between the MPI point and the actual location of the early 
correlator 114 is then found. The difference between ideal 
distance 106 and actual distance 112 is the multipath track- 
ing error 116. 

[0058] Note that because the DLL correlator spacing is 
constant, distances 106 and 112 can be measured from either 
of the two correlators. They can also be measured from the 
prompt correlator at the inferred peak. However, the dis- 
tances must be measured from an MPI point to a point that 
is multipath dependent, i.e., a point whose location indicates 
the performance of the DLL tracking in the presence of 
multipath. This point is referred to herein as a multipath - 
dependent point. Of course, ignoring receiver finite band- 
width effects, the ideal and corresponding actual MPI points 
are at essentially the same location relative to the LOS peak, 
while the ideal and corresponding actual multipath-depen- 
dent points are not at the same location relative to the LOS 
peak. 

[0059] A method according to the invention for using the 
MPI points to improve line-of-sight code tracking has the 
following assumptions: 

[0060] 1. For an individual signal, the relative distance 
between the main lobe and side lobes remains constant. 
That is, the C/Acode can be delayed as a whole (e.g., 
by the ionosphere), but the chip length remains con- 
stant. 

[0061] 2. The receiver identifies the PRN codes before 
beginning the MPI search, and therefore searches for a 
known MPI point. 



[0062] 3. The only interference sources present are 
thermal noise, multipath, and cross-correlation noise 
due to the presence of other visible GPS satellites. 

[0063] 4. The relative positions of receiver correlators 
are controlled by software. 

[0064] FIG. 10 is a flow diagram of a method 120 of the 
invention for mitigating the effects of multipath on code 
tracking of a line-of-sight signal. Method 120 is performed 
for each receiver channel, i.e., for each PRN code. The 
method has an offline portion 122 and an online portion 124; 
Offline portion 122 computes ideal distance 106 (FIG. 9); 
while online portion 124 computes actual distance 112 based 
on the received signal. While the method is illustrated as 
requiring both portions, offline portion 122 is typically 
performed only once for each receiver channel (or PRN 
code). The value of the ideal distance is stored, for example, 
in a lookup table, and then used by online portion 124 to 
determine and mitigate the effect of multipath on tracking of 
the newly received signal. Thus online portion 124 is 
performed continually as signals are received. 

[0065] Beginning with offline portion 122, in step 126, an 
ideal autocorrelation function is computed for the PRN 
code. For example, FIGS. 7A and 7B show ideal autocorr 
relation functions for PRN 8 and 29, respectively. Next, in 
step 128, the MPI point or points are located on the 
autocorrelation function. In the case described above, MPI 
points are located at the late end of zero-slope segments 
(plateaus) of the function that extend over three code chips. 
In general, many MPI points exist within a single autocor- 
relation function, and only the most preferred MPI points are 
selected. For example, MPI points are preferably selected 
from MPI regions with the longest plateaus. Such MPI 
points are valid for a much longer delay in multipath. As 
mentioned above, long-delay multipath with an amplitude 
above the noise floor is unlikely. For each selected MPI 
point, the ideal distance between the MPI point and a 
selected multipath-dependent point is measured (step 130) 
and stored for future reference. The multipath-dependent 
point is typically one of the DLL correlators. As shown in 
FIG. 9, the location of the early or late DLL correlator is 
one-half of the correlator spacing from the ideal autocorre- 
lation peak. 

[0066] Online portion 124 contains analogous steps for the 
received signal. In step 132, the actual correlation function 
between the received signal and the PRN code is computed; 
Note that the PRN code is identical to the PRN code used to 
generate the ideal autocorrelation function in step 126. In 
step 134, one or more multipath-invariant points is identified 
in the computed actual correlation function. The located 
point or points correspond to the points located in the ideal 
autocorrelation function in step 128. Step 134 can be imple- 
mented in many ways; a preferred implementation of step 
134 is described in detail below. Next, in step 136, the actual 
distance between the MPI point and the actual multiparty 
dependent point is measured. The actual multipath-depen- 
dent point must correspond to the multipath-dependent point 
used in step 130. In step 138, the difference between the 
actual distance and the ideal distance is computed. This 
difference, the multipath tracking error, is used in step 140 
to mitigate the effects of multipath on code tracking of the 
line of sight signal by the DLL. Step 140 can also be 
implemented in a variety of ways. 
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[0067] A simple implementation of step 140 is to subtract 
the multipath tracking error from the results of the DLL, i.e., 
from the computed pseudorange to the particular satellite. 

[0068] A more complicated implementation uses the com- 
puted tracking error to control the positions of the DLL 
correlators. Use of the MPI-based multipath corrections to 
adjust the DCO of the DLL may reduce the maximum 
tracking excursions of the DLL due to multipath or thermal 
noise. Accordingly, it may help further reduce the maximum 
expected tracking errors. For example, under conditions of 
significant interference, the DLL correlators can be trans- 
lated so far from the LOS peak that they are unlikely to be 
able to return to the correct position to track the peak. 
Knowing approximate bounds on the distance between the 
MPI point and the LOS correlation peak allows for main- 
taining the correlators within a reasonable distance from the 
peak. 

[0069] Note that the present invention can detect and 
eliminate multipath error introduced by arbitrarily small 
multipath delays, even below 25 m. This is because the 
ability to locate a multipath-invariant point is independent of 
multipath parameters such as delay. At such small relative 
multipath delays, most existing wideband techniques have 
difficulty separating the multipath and line-of-sight main 
lobes and, therefore, perform only as well as narrowband 
techniques. 

[0070] One of the most challenging aspects of implement- 
ing the present invention accurately is locating the multi- 
path -invariant point reproducibly. In general, the signal 
strength is so low at the multipath-invariant point that the 
correlation function has a value that is on the order of the 
noise floor, and it is difficult to extract meaningful data from 
this part of the signal. The problem is particularly difficult in 
receivers with narrow precorrelation bandwidths (PCBs). In 
such receivers, the correlation function has much smoother 
slope transitions than in wider bandwidth receivers. For the 
above example, in which the multipath-invariant property is 
the slope of the correlation function, and the multipath- 
invariant point occurs at a slope transition, narrow PCB 
receiver present particular difficulties. However, because 
few existing multipath mitigation techniques are effective in 
these narrowband receivers, application of the present 
invention to these conditions is particularly desirable. 

[0071] FIG. 11 illustrates the difficulty caused by narrow 
PCB receivers in the method of the invention. A correlation 
function 150 generated by a receiver with infinite PCB has 
very sharp slope transitions, and the location of a multipath- 
invariant point 152 is easily distinguished. A correlation 
function 154 is generated by a receiver with a narrow PCB. 
Function 154 is much smoother than function 150. The 
multipath-invariant point in function 154 that corresponds to 
point 152 is much more difficult to locate. If the point is 
assumed to occur at a value of the correlation function that 
is equivalent to the value at point 152, i.e., at point 156, then 
an error of length 158 is introduced. It is important that the 
criterion for selecting the MPI point in the received corre- 
lation function minimizes error 158. 

[0072] FIG. 12 is a block diagram of a system 160 for 
implementing multipath mitigation in one receiver channel 
according to the present invention. Module 162 implements 
the correlation function computation using the code 
extracted from a modulated carrier signal 164 and the ideal 



code signal for the particular receiver channel. Module 166 
contains the delay-lock loop discriminator that tracks the 
received correlation function. As illustrated, the location of 
the DLL correlators is shifted until the difference between 
early and late samples of the correlation function is zero. The 
location of the tracked peak is passed to navigator module 
168 that computes the pseudorange for this single channel 
and the position solution using pseudoranges from all chan- 
nels. Module 170, which implements multipath mitigation 
according to the present invention, is known as the Tracking 
Error Compensator (TrEC). TrEC 170 has as input the code 
correlation functions 162 and the pseudorange solutions; 
which, combined with the ideal distance from the lookup 
table, can be used to bound the MPI search, as explained 
below. The tracking error value obtained from TrEC 170 is 
subtracted from the pseudorange solution to obtain corrected 
pseudoranges. The tracking error is also passed to navigation 
module 168, which computes the position solution and may 
also perform carrier smoothing of the pseudoranges. Carrier 
signals 164 may be used for carrier smoothing of the 
corrected pseudorange values. Carrier smoothing is a stan- 
dard prior art technique for stabilizing computed position 
and pseudorange values and will not be discussed herein. 
Using the TrEC-corrected pseudoranges results in an MPL 
corrected, carrier-smoothed position solution 165. 

[0073] The specific implementation of TrEC 170 depends 
upon the type of receiver used. For example, TrEC can be 
implemented in software in an existing low-end narrow PCB 
receiver. Obtaining accurate results is most difficult in such 
an implementation, because the correlation function does 
not have sharp peaks, and because limited sampling of the 
correlation function for locating the MPI point is permitted; 
However, this embodiment is relatively inexpensive to 
implement and may not require any additional hardware 
than is provided by conventional receivers. FIG. 13 is a 
block diagram illustrating the elements of TrEC module 170 
of FIG. 12, according to a software-only embodiment of the 
method. TrEC module 170 computes the tracking error from 
the correlation function of the signal, given the results from 
the DLL. 

[0074] Module 172 is an MPI region discriminator, which 
is described in more detail below. In the embodiment 
described below, the region of the correlation function likely 
to contain the MPI point is sampled at two points, and the 
difference between the two points is formed. 

[0075] At each sampling, the value is passed to a filter 
H 174 which time averages the samples. The correlation 
function sampling occurs very rapidiy, at the rate at which 
correlation functions are generated for the DLL (e.g., 10 
Hz). The time-averaged values from filter 174 are passed to 
an optimizer 176 that will also be described in more detail 
below. Briefly, optimizer 176 compares the received value to 
a threshold to determine whether the MPI has been correctly 
identified. If not, it determines how to shift the sample 
locations to approach the MPI point. The result of the 
determination is sent to a digitally controlled oscillator 178 
that controls the location of the MPI correlators 180, i.e., the 
location of the MPI region at which samples are taken. 

[0076] In TrEC 170, carrier dynamics 182 are used to 
decouple the multipath-varying code dynamics from the 
sampled points in the multipath-invariant region. The mul- 
tipath-varying code dynamics are manifested in the move- 
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ment of the DLL correlators. If the multipath region sam- 
pling is not decoupled from the code dynamics, then the 
location of the multipath correlators moves with the DLL 
correlators, which decreases the efficiency and accuracy of 
locating the multipath-invariant point. GPS receivers typi- 
cally perform carrier aiding, in which the dynamics of the 
DLL are propagated from one update to the next using the 
difference in accumulated carrier phase measurements 
between two updates. Thermal noise and other error sources 
affect code phase measurements much more than carrier 
phase measurements, and so carrier aiding effectively 
smoothes code tracking errors. In this context, the carrier 
signal is referred to as a low-distortion component of the 
signal, while the code, from which the correlation function 
is generated, is a high-distortion component. In TrEC 170, 
carrier aiding enables the MPI correlator pair to maintain a 
code-phase offset relative to the DLL that is only carrier- 
phase dependent. 

[0077] Note that while the carrier signal is the preferred 
low-distortion component that is used to decouple dynamics 
from the code signal, this function can be performed using 
any knowledge of the signal dynamics that is not subject to 
the errors affecting the code signal. For example, an inertial 
sensor such as an accelerometer, a separate reference signal, 
or a separate information source can be used. 

[0078] Alternatively, if the distortion of the code dynamics 
is predictable, then the location of the multipath sampling 
can be adjusted based on the predicted code dynamics. In 
this case, the prediction is not a function of any measured 
signal dynamics. 

[0079] When optimizer 176 determines that the MPI point 
has been located correctly, it computes the distance between 
the MPI point and some fixed point of the DLL, such as the 
early correlator or the prompt correlator, using the code 
phase measurements 184 obtained from the DLL. The 
tracking error value is passed to a filter 1^186 that smoothes 
the values over time using carrier phase dynamics 182. Filter 
186 is needed because optimizer 176 makes decisions based 
on average values after a number of samples of the discrimi- 
nator function. As a result, the converged estimate of the 
MPI point location contains some error, which can be 
reduced by filtering over time. Filter 186 is a closed loop 
estimation filter that propagates estimates at one time step to 
the estimate at the next time step using the integrated carrier 
phase. There is no latency associated with dynamic decou- 
pling, and so the solutions can be used immediately to 
compensate for DLL tracking errors. Slow updates (e.g., 
well below 1 Hz) from optimizer 176 are blended with fast 
updates of carrier dynamics 182 (e.g., up to 1 kHz). A 
subsequent estimate of the MPI point is averaged with the 
integrated-carrier-phase-propagated previous solution using 
a Hatch filter with a time constant that can be relatively long, 
because the relative distance to the multipath-invariant point 
is constant. The measured tracking error 188 is then used to 
correct the pseudoranges. 

[0080] In this embodiment, MPI region discriminator 172 
is implemented using a single correlator pair; that is, the 
correlation function is sampled at two points that are typi- 
cally a fixed distance apart. The goal of the single-correlator 
implementation of the MPI discriminator is to find a point at 
a particular threshold value of the correlation function. This 
threshold value defines the MPI point. FIG. 14A indicates a 



multipath-invariant region 190 of a correlation function 
containing line-of-sight 191 and multipath 193 contribu- 
tions. In this and remaining figures, x is the time relative to 
the line-of-sight peak, and T c is the chip length in time units. 
As described above, the preferred multipath-invariant point 
is located at the slope discontinuity 192. Of the points in 
multipath-invariant region 190, the point at location 192 is 
easiest to distinguish and remains multipath invariant over 
the longest multipath delay. In a narrow PCB receiver, the 
correlation function does not have sharp peaks but looks 
more like the curve of FIG. 14B. In this case, the equivalent 
multipath-invariant point 200 is located in a region of 
non-zero slope. In general, even if the correlation function 
has sharp corners, the multipath-invariant point is easiest to 
locate if it is considered to be a small distance up the side of 
the main lobe. This position is referred to as being adjacent 
to a plateau of the correlation function. Referring again to 
FIG. 14A, the point must be in a region 194 between the 
slope transition at location 192 and the beginning of the 
multipath peak at a location 196. Assuming that the multi- 
path delay always has a value of at least 6, region 194 is 
guaranteed to be multipath invariant, while any point to the 
right of location 196 is not multipath invariant. To minimize 
the assumed value of 8, the location at which the multipart^ 
invariant point is assumed to exist should be within region 
194 and close to location 192. A value e of the function 
within region 194 is selected as the threshold value used to 
locate the MPI point on the received correlation function. 

[0081] Note that it is the slope and not the amplitude that 
is multipath invariant. If a threshold value is to be used to 
identify the multipath-invariant point, it is important to 
remove the effects of multipath on the correlation function 
amplitude, so that only the slope is effectively being mea- 
sured. This is done by forming a difference between the 
correlation function value at two locations, e.g., an early 
correlator 202 and a late correlator 204, as illustrated in FIG. 
14B. 

[0082] The difference between the sampled value at cor- 
relator locations 202 and 204 is referred to as a discriminator 
function, which is then compared with the predefined thresh-- 
old value. Assuming that the multipath-invariant region 
contains no line-of-sight side lobes, the value of the total 
correlation function sjftjs^l) at point 202 (early) is 

[0083] where b 0 is a noise-level constant value and p is the 
multipath contribution, which is a function of multipath 
amplitude, phase, and relative frequency. The value of the 
total correlation function at point 204 (laie) is 

[0084] where R(t-tJ is the LOS contribution from the 
main lobe at a point x c . Note that points 202 and 204 are 
within the multipath-invariant region. They may have a 
multipath component, but, by definition, the multipath comr 
ponent is constant in this region, i.e., independent of mul- 
tipath parameters. Thus the p terms at the two locations are 
identical, and the difference between the correlation function 
at the two points is simply given by R^-tJ; the constant 
value b 0 and the multipath contributions cancel. The disr 
criminator function therefore indicates how far up the line- 
of-sight main lobe late correlator 204 is. This quantity is 
exactly what is required to locate the multipath-invariant 
point at the chosen threshold value. To implement this 
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embodiment, a threshold value is selected for the difference 
between late 204 and early 202 correlators, and the correla- 
tor position is adjusted until the discriminator value is equal 
to the threshold value, within an allowed tolerance. 

[0085] In FIG. 13, filter 174 time averages discriminator 
values. Optimizer 176 determines, based on the difference 
between the discriminator value and the predetermined 
threshold e, in which direction and by how much to move the 
correlator pair so that the late correlator is positioned at the 
multipath-invariant point. Preferably, optimizer 176 is 
implemented using a Fibonacci search. A Fibonacci search 
is an algorithm for locating the minimum (or maximum) of 
a function within a predefined interval of length I. Similar to 
a binary search, a Fibonacci search continually narrows the 
interval until the desired point is located to within a given 
tolerance. In subsequent iterations, the interval bounds and 
test points within the interval are adjusted according to a 
Fibonacci sequence, in which Fj"^-! and F^F^+F^. 
For a discrete interval, the number of iterations is set by the 
number of points within the original interval. For a continu- 
ous interval, the number of iterations is set by the desired 
precision, i.e., the size of the final interval. The final interval 
size is given by I/F n , where n is the number of iterations. A 
typical Fibonacci search to minimize a function f(x) within 
an interval [b^bj proceeds as follows: 

[0086] Initialization: 
[0087] j th Iteration: 

if Ax^m^- 

b u -x 2 

X 2 «'X 1 

^l=(Vfn-j-A-j)Ml) 

if/(*iM*2): 

bj-X! 

[0088] The iterations are repeated through the j=n-l itera- 
tion. The solution is located within the final interval, i.e., 
between b 1 and b u , 

[0089] A Fibonacci search has a number of properties that 
make it suitable for locating the multipath-invariant point. 
Because it uses the Fibonacci ratio to reduce the interval, 
rather than halving the interval as a binary search does, it is 
a very efficient, non-gradient search for a large initial 
uncertainty. It is also well suited for constrained (bounded) 
uncertainty intervals and functions that are costly to evaluate 
or observe. Each additional iteration of sampling the corre- 
lation function increases the initialization time to locate the 
multipath-invariant point, and so a method that rapidly 
zeroes in on the point is highly desirable. 

[0090] In this embodiment of the method, a standard 
Fibonacci search (FS) is adapted to the search for the 
multipath-invariant point. Rather than searching for a mini- 
mum of a function in a fixed interval, the search is for a point 
yielding a particular value, the threshold e, of the discrimi- 
nator function within a fixed interval. The discriminator 
function has a value of Late -Early, where Late and Early are 
the correlation function samples at points x mpi and (x mpi - 



d early ) respectively, and where T mpi is the offset from the 
main peak at which the late correlator is sampled and d early 
is the time difference between the two correlator locations: 

[0091] The other significant difference between the MPI 
point search and a standard Fibonacci search is that the 
correlation function is dynamic. The dynamics are removed 
as described above. Because evaluations of the correlation 
function are averaged over a period of time, the correlator 
pair maintains a fixed code offset during averaging. 

[0092] The FS is initialized with an initial guess at x mpi ; 
Initially, the location of the MPI point relative to this guess 
is unknown. All subsequent FS updates (i.e., decisions and 
movements) are made relative to this position. Again, all 
dynamics are removed at each step using carrier phase 
propagation. Each decision narrows the interval until a final 
estimate of T mpi is made. This final estimate is then stored 
and propagated in time using the carrier phase. Indepen T 
dently, the FS is repeated and a second guess of the MPI 
point location (x mpi ) is made. This guess is averaged (filtered 
using Hp) with the first to produce an improved estimate. 
The improved estimate is then propagated forward and the 
process is repeated. The result of the search is a value of the 
late correlator offset from the main line of sight peak, 
that yields a discriminator value of e, the chosen threshold; 
This point is defined to be the multipath-invariant point. : 

[0093] To find the multipath-invariant point, the steps are 
as follows: 

[0094] Initialization: 
[0095] j* Iteration: 

if Late(x mpi )-Early(T mpi -rf (;illy )^€: (multipath-invari- 
ant point to the right) 

if Late(T a)pi )-Eailv(T mpr d eaily )>e: (multipath-invari- 
ant point to the left) 

[0096] The iterations are repeated through the j=n-l itera- 
tion. The multipath-invariant point is located within the final 
interval, i.e., between b u and b a . For example, for a narrow- 
band receiver with a DLL correlator spacing of 0.5 chips; 
I«±60 meters yields F n =144. This implies 11 evaluations and 
final interval size ^1 meter. The first Fibonacci search has 
a relatively large initial interval, but subsequent iterations 
will preferably have smaller initial intervals. 

[0097] Initial values of lower and upper bounds b a and b„ 
of the Fibonacci search interval, shown in FIG. 14C, are 
determined from the current position of the delay lock loop 
prompt correlator 212 (i.e., the midpoint between the early 
and late correlators), the range of expected multipath-in- 
duced tracking error E, and the ideal distance 214 between 
the multipath-invariant point 216 and DLL correlator 210 of 
the ideal autocorrelation function. Suitable values of the 
error E are obtained from known envelopes (FIG. 6) bound- 
ing the multipath tracking error of a standard DLL. For 
example, for a 0.5-chip spacing DLL and a signal-to^ 
multipath ratio of 3 dB, multipath tracking error has a 
maximum range of approximately -50 m to +50 m, equiva- 
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lent to 0.17 chips earlier and later than the ideal distance. 
The initial search interval bounds are computed as: 

^p«Hnpr- rf id«d-£ and 

[0098] where d^i is ideal distance 214, +E is the tracking 
error, and x prompt is the actual location of the DLL prompt 
correlator 212, i.e., the point from which the ideal distance 
is measured. 

[0099] In an alternative embodiment, a ratio of two corr- 
elator values is used instead of a difference between two 
values, i.e., instead of the discriminator function. In this 
embodiment, the ratio of two values is compared with a 
predetermined threshold value that indicates the location of 
the multipath-invariant point. The threshold value is typi- 
cally different from a threshold value selected based on 
discriminator values; however, all other implementation 
details are the same. The Fibonacci search is performed 
using this ratio rather than the discriminator function. The 
reasoning behind the ratio is similar to the reasoning behind 
the discriminator function. Consider that each correlator 
function sample point is in the form Ax, where A is a 
multipath-dependent scale factor and x is a line-of-sight 
contribution. Since samples are taken in the multipath- 
invariant region, A is constant for all the samples. A ratio of 
two samples is therefore multipath independent and indi- 
cates the location of the point relative to the line-of-sight 
function. 

[0100] Referring again to FIG. 13, the final result of the 
Fibonacci search is passed from optimizer 176 to filter 186. 
Because the Fibonacci search converges only after a pre- 
scribed number of iterations, the rate at which a solution is 
passed from optimizer 176 to filter 186 is much slower than 
the rate at which correlation function samples are obtained 
in discriminator module 172. 

[0101] One drawback of this embodiment is the relatively 
long time required to initialize the MPI location process. 
Because each iteration of the Fibonacci search samples only 
one location of the correlation function, the MPI point and 
resulting tracking error value are only sufficiently accurate 
after a period of minutes. The Fibonacci search using the 
difference between values of the correlation function at two 
different points (late minus early) is the preferred embodi- 
ment for implementing the present invention on a low-end 
receiver in which only a single additional correlator pair is 
available. Other methods are preferable if additional hard- 
ware can be added to the receiver. In these embodiments, the 
correlation function is preferably sampled at multiple points 
to locate one or more multipath-invariant points. Additional 
hardware embodiments have the advantages of lower ini- 
tialization time and/or higher accuracy. 

[0102] FIG. 15 illustrates a curve fit multiple correlator 
embodiment. As with the previous embodiment, it is to be 
understood that multiple correlator pairs represent either 
additional hardware or, in purely digital receivers, multiple 
sampling near the multipath-invariant point. In this embodi- 
ment, the correlation function is sampled at multiple early 
and late positions. If the correlators are separated by a fixed 
distance, then each of the early correlators is within a region 
that is unaffected by the main lobe of the line-of-sight 
correlation function. Alternatively, if samples can be 
acquired at arbitrary locations, then one early sample is 



taken and subtracted from each of the late samples to arrive 
at values of the discriminator function at each late sample 
time. FIG. 15 illustrates four different sample points, t 1 
through x 4 , although any number of sample points can be 
used. A curve fit is generated from either the sample points 
or the discriminator function. The point at which the curve 
fit is at the threshold value e is considered to be the 
multipath-invariant point. If a sufficient number of correlar 
tors are provided, then no search is required, and the MPI 
point is located based on a single curve fit. In this case, 
continuous updates to the MPI point estimate can be made 
and smoothed over time. Note that this embodiment effec- 
tively acquires the same amount of information in one 
sampling time as the single-correlator pair embodiment 
acquires in four sampling times. That is, if each of points x^ 
through x 4 requires sampling at a different time in the 
single-correlator embodiment, such an implementation will 
take at least 4 times as long to estimate the MPI point. 

[0103] FIG. 16 illustrates an alternative multiple correlar 
tor embodiment. In this embodiment, the Fibonacci search 
algorithm is implemented for a single correlator pair located 
at each of a plurality of multipath-invariant points. For 
example, a first multipath-invariant point is located at an 
oflfeet of x v The difference between an actual distance 220 
and an ideal distance 222 is a tracking error 224. A second 
multipath-invariant point is located at an offset of x 2 . 

[0104] The difference between a second actual distance 
226 and a second ideal distance 228 is a tracking error 230; 
A third multipath-invariant point is located at an offset of x 3 ; 
The difference between a third actual distance 232 and a 
third ideal distance 234 is a tracking error 236. In theory, as 
illustrated in FIG. 16, all errors 224, 230, and 236 have the 
same value. In principal, however, if one of the multipath 
invariant points is not located correctly, there is a small 
variation in the errors. Individual tracking errors 224, 230, 
and 236 are averaged to obtain an overall tracking error that 
is used to correct the pseudorange. One advantage of this 
embodiment is an improvement in the accuracy of the 
average tracking error compared with individual tracking 
errors derived from a single multipath-invariant point; 
Assuming the optimization errors are randomly distributed 
about the true MPI point location, they can be reduced by a 
factor of the number of MPI correlator pairs per channel. A 
second advantage is that the individual MPI estimates can be 
improved by using the code phase differences between MPI 
correlators. These differences are affected by the TrEG 
implementation bias errors (discussed below) only. Com- 
paring all code phase offsets allows outliers to be discovered 
and removed. This embodiment is implemented either by 
using additional hardware or by sampling the correlation 
function at locations near each multipath-invariant point. ; 

[0105] The two multiple-correlator embodiments can also 
be combined by using multiple correlators to sample each of 
a number of different MPI points. 

[0106] The multipath mitigation method of the present 
invention has been implemented using signal generators 
with known multipath delays and using actual satellites. In 
all cases, the mean tracking error is reduced significantly! 
FIG. 17 shows theoretical multipath performance of Xh6 
present invention, as well as existing narrowband and wide- 
band techniques, implemented on a 2 MHz PCB receiver 
with a signal-to-multipath ratio of 3 dB. Note that unlike the 
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performance envelopes for existing techniques, tracking 
error of the present invention displays a constant bias. It is 
constant because the multipath error that occurs with the 
present invention is not a function of multipath parameters. 
The bias results from the finite PCB of the receiver and the 
selection (location) of the desired MPI point. TrEC bias 
errors, which appear as relatively small variations about this 
constant bias, result from a combination of five factors: 

[0107] 1. Optimization tolerance. The Fibonacci search 
converges on the MPI point within a predetermined 
tolerance. Without noise, the error can be made arbi- 
trarily small; with noise, however, the tolerance must 
be increased to provide reasonable initialization times. 

[0108] 2. Thermal noise. Low signal power in the MPI 
region makes thermal noise significant. The biases are 
higher for low signal power, low-elevation satellites 
and with large search intervals. Long integration times 
decrease the effect of thermal noise. 

[0109] 3. Very short delay multipath. 

[0110] 4. Amplitude variations. 

[0111] 5. Ionospheric divergence rate. Low-elevation 
satellites cause the largest errors, which are magnified 
by long integration times and long-term smoothing. 

[0112] One way to remove these biases is by additional 
hardware (e.g., wider receiver precorrelation bandwidth or 
more correlators). In addition, because the tracking error 
bias is multipath independent, its value can be determined 
once and then used to correct subsequent solutions. Various 
methods can be used to determine the bias. For example, the 
frequency band limiting effects of a particular receiver can 
be modeled using a PCB filter model, and the resultant bias 
in the tracking error from these effects can be determined. 
Alternatively, the bias can be found directly by using a GPS 
signal generator or pseudolite. Provided the receiver's con- 
ventional DLL is not corrupted by multipath, the bias can be 
computed as simply the mean difference between the ideal, 
infinite bandwidth distance and the measured distance to the 
MPI point. The computed bias is then added to subsequent 
tracking error measurements. 

[0113] Note that changes in the satellite signal power 
change the amplitude of the correlation function and the 
location at which the selected threshold value occurs. This 
causes the filter bias to be different for each satellite. 
Because the received signal power is a function of the 
satellite elevation angle, the correlation peak is also eleva- 
tion dependent. This effect can be eliminated by normalizing 
the correlation function. Preferably, the Fibonacci search is 
performed using a normalized correlation function. A nor- 
malized MPI discriminator D'(e) is given by: 



[0114] where 2(I lfinpi _ Ic mpi ) is the integrated mean differ- 
ence of the late and early in-phase samples in the MPI 



region, V^XQiw+Qe^i) is the accumulated late and early 
quadrature samples, I P is the in-phase prompt measurement 
obtained at the DLL, and Q P is the quadrature prompt 
measurement obtained at the DLL. 

[0115] If the signals are accordingly normalized, not only 
is the tracking error bias independent of multipath paramr 
eters such as amplitude, phase, and delay, it may also be 
independent of PRN code. That is, the bias is constant across 
all receiver channels for all LOS signal variations. This 
provides effective normalization with respect to signal 
power variations, provided the multipath and the LOS signal 
attenuation effects can be modeled as simple amplitude scale 
factors. This also assumes the accumulated quadrature 
samples, ZQ p and Vi2(Q limp£ +Q e(inpi ), are nonrandom and 
possess nonzero mean. 

[0116] With a constant PCB filter (and MPI point location) 
bias, even though the bias affects the individual pseudorange 
measurements, its effect can be removed from the final 
position solution. If the bias is the same across channels, it 
becomes part of the clock bias in the final GPS position 
solution. 

[0117] The filter bias can also be reduced through calibra- 
tion by removing the elevation angle dependence of the 
received signal power. Beam shaping of the transmitted GPS 
signal has a relatively minor effect, causing the line-of-sight 
signal power to vary by only 2-3 dB as a function of 
elevation angle. The dominant factor causing variation is the 
antenna gain pattern, which varies widely for different 
antenna types. Once known, the elevation angle dependence 
is removed by adjusting the MPI tracking threshold as a 
function of satellite elevation angle and, if the pattern is not 
symmetric, azimuth. Calibration is simply experimental 
determination of a scale factor a(8 elev ) such that a(9 e i ev )=A* 
for k ranging over all satellites, where A is the line-of-sight 
signal amplitude. The calibration-normalized discriminator 
function is the discriminator function scaled by a(G elev ), with 
the threshold scaled similarly. 

[0118] Cross-channel smoothing (CCS) is a technique for 
reducing the TrEC bias errors and can be used to reduce the 
individual pseudorange biases. In CSS, the tracking errors 
are averaged across all receiver channels. If the nominal 
TrEC correction for the VP pseudorange at time tj is defined 
as Wk(tj), a modified, cross-channel smoothed TrEC correc- 
tion, *P xc (tj), may then be defined as 

[0119] where 



i K 



[0120] is the average of all K original TrEC corrections at 
time tj. ^ k (tj) represents the time-averaged (Hatch- filtered) 
estimate of the k* TrEC corrections at time tj. Cross-channel 
smoothing effectively reduces the maximum TrEC bias 
errors by a factor of the number of TrEC-enabled channels 
on a GPS receiver. As a consequence of this, however, the 
TrEC corrections become multipath varying, and the maxi^ 
mum expected multipath error increases. This trade-off; 
however, is acceptable in many instances where small mean 
errors are most desirable. Cross-channel smoothing may not 
be necessary if more MPI correlators and/or a wider PCB is 
available. 
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[0121] While the invention has been described with 
respect to correlation function slope as the multipath-invari- 
ant property, other multipath-invariant properties may exist. 
Appropriate changes to the algorithm are required to locate 
multipath-invariant points that are characterized by different 
properties. It will be apparent to one of average skill in the 
art how to implement the necessary changes. 

[0122] It is believed that the present invention provides 
one of the few multipath mitigation techniques that are 
effective in narrow precorrelation bandwidth receivers. Is it 
therefore possible that the invention will be implemented on 
a widespread basis. If so, it will be helpful to use code 
sequences that are optimal for locating multipath-invariant 
points. Optimal codes (e.g., m-sequenccs) have autocorre- 
lation functions with long multipath-invariant regions adja- 
cent to the main lobe. 

[0123] While the present invention has been described 
primarily for use in mitigating multipath, it can be used for 
detecting and eliminating tracking error caused by any 
source. An analogous algorithm is used, with multipath- 
invariant points replaced by points that are distortion invari- 
ant, i.e., that have a property that is independent of param- 
eters of the signal distortion. One important application is 
integrity monitoring, i.e., monitoring the integrity of the 
transmitted signals to detect signal anomalies. In high- 
precision GPS applications such as landing aircraft, a very 
high integrity of the transmitted signal is required. Satellite 
signal anomalies result from a failure of the signal gener- 
ating hardware on one of the GPS space vehicles (SVs). For 
example, a failure of SV19 occurred in October 1993 and 
caused pseudorange errors on the order of 3 to 8 meters, 
which is much too large for aircraft conducting a precision 
approach. The signal integrity is monitored independently of 
its direct use for navigation, and detected errors in signal 
transmission are removed from the position estimate. 

[0124] Signal transmission failures distort the correlation 
peak in a manner similar to the effect of multipath. Rather 
than a multipath-invariant point, however, points that are 
insensitive to satellite failure are first located on the ideal 
autocorrelation function. During operation, these points are 
located in the received correlation function using methods 
identical to those described above for locating multipath- 
invariant points. The difference of the distances between the 
ideal and actual failure -insensitive points and the DLL loop 
is used as either an indicator of signal failure or as com- 
pensation for tracking error caused by signal failure. These 
points may or may not be coincident with MPI points. 

[~0 1 ^ ^1 A nntK^r 1 m -r\r\ rfint qnnli'r>nt(/Mi r\ f inimnti^p Ir 

|_VJ.*»fc>J 1 lilUlUVl ILXA^JKJ L I CI 14 L 1 1VU Ki\J 11 Ul 111V UlYWiiliOLl lO 

mitigating tracking error caused by noise in the received 
code. Multipath-invariant points are also insensitive to code 
noise. In fact, experiments performed by the present inven- 
tor indicate that the invention is capable of eliminating both 
multipath tracking error and code noise tracking error. In 
practice, it is not necessary to separate the two effects, as the 
effect on the correlation peak is qualitatively the same. For 
systems in which significant multipath is not present, the 
present invention can be implemented for code noise track- 
ing error reduction alone. Implementation is identical to that 
described above. The chief advantage comes from the 
dynamic decoupling enabled by the (low-noise) carrier 
phase tracking. Tracking the MPI point corrects for any code 
tracking biases, while the relatively stable carrier dynamics 
effectively cancel the random noise errors. 



[0126] In applications in which multipath is not the pri- 
mary source of code tracking error, the distortion-invariant 
point may be located at the late (i.e., right) side of a peak or 
side lobe of the correlation function. While multipath dis- 
torts the late side of peaks, other error sources distort the 
early side of the peak. For example, signal transmission 
errors may cause signal anomalies known as evil wave- 
forms. One type of evil waveform, Threat Model A of the 
full ICAO second-order step threat models, is a lead/lag 
model that may cause signal distortions that lead the signal: 
In this case, the distortion-invariant point is at the late side 
of peaks. 

[0127] The present invention may be implemented to use 
either a combination of correlators (e.g., using cross-corre- 
lation properties) or a specially-designed non-replica code 
function to alter the received correlation function. Sampling 
with such codes may change the characteristics of the 
correlation function, thereby providing additional signal 
power for more robust solutions in the presence of noise. 
Such technology is used in other multipath mitigation tech- 
niques and can be applied to the present invention. For 
example, the Strobe Correlator is described in L. Garin and 
J. -M. Rousseau, "Enhanced Strobe Correlator Multipath 
Mitigation for Code and Carrier, "Proceedings of the ltf* 
International Technical Meeting of the Satellite Division of 
the Institute of Navigation, Vol. 1, pp. 559-568, 1997. 

[0128] Improved (loop-filter) estimation algorithms may 
assist in achieving faster, more robust convergence to the 
MPI point. A traditional Hatch filter was assumed for the 
TrEC implementation described above; however, using 
SNR, elevation angle, and other signal observables in ah 
adaptive estimator or Kalman filter may provide even more 
effective TrEC corrections. The model equations for these 
additions are available in the prior art and are not described 
herein. In the context of the present invention, these error 
estimates are applied in combination with the TrEC estiT 
mates, which may be particularly useful when hardware is 
limited and when CCS is necessary. It is to be understood 
that the combination of multipath-invariant point location 
according to the present invention with prior art multipath 
estimation techniques is within the scope of the present 
invention. 

[0129] For GPS signals, the ionosphere typically causes 
the code signal to diverge from the carrier. This limits the 
time that carrier phase propagation (or dynamic decoupling) 
can be used in single-frequency receivers. Dual-frequency 
measurements can be used, however, to estimate this rate of 
divergence. In GrPS, this rate m formation may be obtained 
from dual-frequency receivers or from WAAS (Wide Area 
Augmentation System) receivers. If this rate is available, it 
may be used directly to correct the carrier propagation rate 
and accordingly permit longer sample averaging times and/ 
or longer propagation times. 

[0130] It will be clear to one skilled in the art that the 
above embodiments may be altered in many ways without 
departing from the scope of the invention. Accordingly, the 
scope of the invention should be determined by the follow-i 
ing claims and their legal equivalents. 

What is claimed is: 

1. A method for mitigating the effect of multipath signals 
on envelope tracking of a primary signal by a spread- 
spectrum receiver, wherein said primary signal comprises a 
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carrier signal modulated by an envelope, and wherein a 
received signal comprises said multipath signals and said 
primary signal, said method comprising: 

a) locating an ideal multipath-invariant point in an ideal 
comparison function, thereby determining an ideal dis- 
tance between said ideal multipath-invariant point and 
an ideal multipath-dependent point, wherein said ideal 
comparison function is computed from an ideal signal 
that does not comprise multipath signals; 

b) locating a corresponding multipath-invariant point in 
an actual comparison function, thereby determining an 
actual distance between said corresponding multipath- 
invariant point and a corresponding multipath-depen- 
dent point, wherein said actual comparison function is 
computed from said received signal; 

c) comparing said ideal distance and said actual distance; 
and 

d) based on said comparison, mitigating the effect of said 
multipath signals on said envelope tracking of said 
primary signal. 

2. The method of claim 1 wherein step (c) comprises 
subtracting said ideal distance from said actual distance to 
obtain a tracking error. 

3. The method of claim 2 wherein step (d) comprises 
adjusting a computed pseudorange by said tracking error, 
wherein said pseudorange is computed in dependence on 
said envelope tracking and represents a distance between 
said receiver and a signal transmitter. 

4. The method of claim 2, further comprising removing a 
bias in said tracking error. 

5. The method of claim 2 wherein said receiver further 
tracks a plurality of additional received signals, and wherein 
said method further comprises: 

performing steps (a), (b), and (c) for each additional 
received signal to obtain a plurality of additional track- 
ing errors; and 

averaging said tracking error and said plurality of addi- 
tional tracking errors to obtain an average tracking 
error. 

6. The method of claim 1 wherein step (d) comprises 
controlling a signal tracking loop that tracks said received 
signal. 

7. The method of claim 1 wherein said ideal multipath- 
invariant point is in a plateau of said ideal comparison 
function. 

8. The method of claim 7 wherein said ideal multipath- 
invariant point is at an edge of said plateau. 

9. The method of claim 1 wherein said ideal multipath- 
invariant point is adjacent to a plateau of said ideal com- 
parison function. 

10. The method of claim 1 wherein said corresponding 
multipath-dependent point is a point sampled by a signal 
tracking loop that tracks said received signal. 

11. The method of claim 1 wherein step (b) comprises: 

bounding an interval containing said corresponding mul- 
tipath-invariant point, wherein bounds of said interval 
are selected in dependence on theoretical tracking 
errors and on said ideal distance; and 

narrowing said interval until said corresponding multi- 
path-invariant point is located. 



12. The method of claim 11 wherein narrowing said 
interval comprises performing a modified Fibonacci search. 

13. The method of claim 1 wherein step (b) comprises 
sampling said actual comparison function at at least one 
point to obtain at least a first sample. 

14. The method of claim 13 wherein sampling said actual 
comparison function comprises eliminating the effect of 
dynamics of a signal tracking loop on said at least one point 
using a low-distortion component of said received signal. 

15. The method of claim 14 wherein said low-distortion 
component of said received signal is said carrier signal. 

16. The method of claim 14, further comprising estimatr 
ing a rate of divergence of said signal tracking loop from 
said low-distortion component. 

17. The method of claim 16 wherein estimating said rate 
of divergence comprises obtaining dual-frequency measure- 
ments of said received signal. 

18. The method of claim 13 wherein step (b) further 
comprises sampling said actual comparison function to 
obtain at least a second sample. 

19. The method of claim 18 wherein step (b) further 
comprises: 

subtracting said first sample from said second sample to 
obtain a discriminator value; and 

comparing said discriminator value with a predetermined 
threshold value. 

20. The method of claim 19 wherein said first sample is 
not affected by a main lobe of a primary signal comparison 
function. 

21. The method of claim 19, further comprising repeating 
said sampling, subtracting, and comparing steps at different 
locations until said discriminator value is sufficiently close 
to said threshold value. 

22. The method of claim 18 wherein step (b) further 
comprises: 

forming a ratio between said first sample and said second 
sample; and 

comparing said ratio with a predetermined threshold 
value, 

23. The method of claim 18 wherein step (b) further 
comprises estimating a shape of said actual comparison 
function from at least said first sample and said second 
sample, and locating said corresponding multipath-invariant 
point in dependence on said shape. 

24. The method of claim 1, further comprising: 

locating one or more additional ideal multipath-invariant 
points in said ideal comparison function, thereby deter^ 
mining one or more additional ideal distances between 
said one or more additional ideal multipath-invariant 
points and said ideal multipath-dependent point; and i 

locating one or more additional corresponding multipath- 
invariant points in said actual comparison function^ 
thereby determining one or more additional actual 
distances between said one or more additional corre4 
sponding multipath-invariant points and said corre- 
sponding multipath-dependent point. 

25. The method of claim 24, further comprising compar- 
ing said one or more additional ideal distances and said one 
or more additional actual distances. 
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26. The method of claim 24, further comprising compar- 
ing said actual distance and said one or more additional 
actual distances. 

27. The method of claim 1, further comprising repeating 
step (b) at a later time to obtain a plurality of additional 
actual distances, and averaging said actual distance and said 
plurality of additional actual distances to obtain an average 
actual distance. 

28. The method of claim 1 wherein said actual compari- 
son function is further computed from a function that is not 
a replica of said received signal. 

29. The method of claim 1, further comprising calibrating 
said actual comparison function in dependence on an eleva- 
tion angle of a signal transmitter. 

30. The method of claim 1, further comprising normaliz- 
ing said actual comparison function. 

31. The method of claim 1 wherein said receiver is a 
Global Positioning System (GPS) receiver. 

32. A method for measuring a tracking error generated by 
a spread-spectrum receiver tracking a received signal com- 
prising an envelope modulating a carrier signal, said method 
comprising: 

a) locating an ideal distortion-invariant point in an ideal 
comparison function, thereby determining an ideal dis- 
tance between said ideal distortion-invariant point and 
an ideal distortion-dependent point, wherein said ideal 
comparison function is computed from an ideal signal 
that is not distorted; 

b) locating a corresponding distortion-invariant point in 
an actual comparison function, thereby determining an 
actual distance between said corresponding distortion- 
invariant point and a corresponding distortion-depen- 
dent point, wherein said actual comparison function is 
computed from said received signal, wherein said 
received signal is at least partially distorted; and 

c) subtracting said ideal distance from said actual distance 
to obtain said tracking error. 

33. The method of claim 32, further comprising using said 
tracking error to improve tracking performance of said 
receiver. 

34. The method of claim 32, wherein said tracking error 
is caused by a failure in transmission of said received signal, 



and wherein said distortion-invariant point has a property 
that is unaffected by said failure. 

35. The method of claim 32, wherein said tracking error 
is caused by noise in said received signal, and wherein said 
distortion-invariant point has a property that is unaffected by 
said noise. 

36. The method of claim 32, further comprising using said 
tracking error to improve tracking of a different received 
signal by said receiver. 

37. The method of claim 32, further comprising removing 
a bias in said tracking error. 

38. The method of claim 32, further comprising calibrat- 
ing said actual comparison function in dependence on an 
elevation angle of a signal transmitter. 

39. The method of claim 32, further comprising normal- 
izing said actual comparison function. 

40. The method of claim 32 wherein step (b) further 
comprises eliminating the effect of dynamics of a signal 
tracking loop on locating said distortion-invariant point. 

41. The method of claim 40 wherein said effect is elimi- 
nated using a low-distortion component of said received 
signal. 

42. The method of claim 41 wherein said low-distortion 
component of said received signal is said carrier signal. 

43. The method of claim 41, further comprising estimate 
ing a rate of divergence of said signal tracking loop from 
said low-distortion component. 

44. The method of claim 43 wherein estimating said rate 
of divergence comprises obtaining dual-frequency measure^ 
ments of said received signal. 

45. The method of claim 32 wherein said receiver further 
tracks a plurality of additional received signals, and wherein 
said method further comprises: 

performing steps (a), (b), and (c) for each additional 
received signal to obtain a plurality of additional track- 
ing errors; and 

averaging said tracking error and said plurality of addi- 
tional tracking errors to obtain an average tracking 
error. 

* * * * * 



07/22/2004, EAST Version: 1.4.1 



